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Hydrogen dynamics and light-induced structural changes in hydrogenated amorphous 

silicon 
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(Dated: February 6, 2008) 

We use accurate first principles methods to study the network dynamics of hydrogenated amor- 
phous silicon, including the motion of hydrogen. In addition to studies of atomic dynamics in the 
electronic ground state, we also adopt a simple procedure to track the H dynamics in light-excited 
states. Consistent with recent experiments and computer simulations, we find that dihydride struc- 
tures are formed for dynamics in the light-excited states, and we give explicit examples of pathways 
to these states. Our simulations appear to be consistent with aspects of the Staebler-Wronski effect, 
such as the light-induced creation of well separated dangling bonds. 

PACS numbers: 61.43.-j, 66.30.-h,71.23.-k,78.20.Bh 
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I. INTRODUCTION 

A variety of experiments and many theoretical stud- 
ies of light-induced structural changes in hydrogenated 
amorphous silicon (a-Si:H) reveal the complexity of pho- 
toresponse in these materials. Light-induced changes in 
photoconductivity and defect format ioni*2^, enhanced 
hydrogen diffusion^*S*L£, and creation of preferential 
proton separation^ are among the phenomena observed 
experimentally in a-Si:H. A special feature of these ma- 
terials is the Jekyll-Hyde behavior of hydrogen as a de- 
fect passivator enabling the practical utilization of the 
material but also as a culprit in light-induced defect 
creationiSiii*i2ii^. A variety of experiments implicate 
H and its dynamics as being important player to the 
SWEil 

There are various models proposed to explain the light- 
induced metastability. Chang et q/i 15 : 16 : 17 suggested that 
the dissociation of a two hydrogen interstitial complex, 
(H|), into separate and more mobile H atoms, caused 
by carriers localized on the H|, is a mechanism for the 
metastable phenomena. In the hydrogen flip model, 
Biswas et al. demonstrated a higher energy metastable 
state is formed when H is flipped to the backside of the 
Si-H bond at a monohydride sitesi^. In the hydrogen 
collision model proposed by Brana&, the recombination- 
induced emission of H from Si-H bond creates mobile H 
and dangling bonds, and these newly created dangling 
bonds become metastable when two mobile H atoms col- 
lide to form a metastable complex containing two Si-H 
bonds. The two-phase model of Zafar and Schiff 19 20 , 
explained thermal stability data, exploited the concept 
of paired hydrogen, and later merged with Branz model 
and invoked dihydride bonding^. Though they have 
different detailed mechanisms for the formation of the 



structures, these models share the same notion that the 
diffusion of H and H pair formation plays a key role 
in light-induced metastability. From ab initio simula- 
tion in a photo-excited state, Fedders, Fu and Drabold 
showed that changes of charge of well-localized defect 
states could induce defect creation2£. Our work can be 
regarded as a natural extension of Ref. |53 to the case of a 
hydrogenated system, and with more accurate techniques 
than those available at the time of the original study. 

Our work has been particularly driven by the experi- 
mental work of Su et al^ who performed proton NMR 
experiments in a-Si:H and found preferential creation of 
H-H distance of 2.3±0.2 A. We have shown that SiH 2 
configurations in the solid state are consistent with these 
observations^. More recently, an NMR study by Bobela 
et aU^, indicated a shorter proton-proton distance in a 
sample of a-Si:H with a somewhat higher defect density. 

In a recent paper—, we briefly reported the hydrogen 
dynamics and light-induced formation of a SiH2 struc- 
ture for a small (71 atom) model. In this article, we 
provide detailed results of simulation on electronic prop- 
erties, vibrational properties and mechanism of hydrogen 
diffusion in both the electronic ground state and light ex- 
cited state in a-Si:H using two different supercell models: 
a larger 223 atom a-Si:H model (Model-I) and a 71 atom 
a-Si:H model (Model-II). We start with small, but topo- 
logically credible models of a-Si:H and use a sufficiently 
accurate method to simulate network dynamics, includ- 
ing the diffusive motion of hydrogen. In a photo-excited 
MD simulation, we find enhanced H motion and the pref- 
erential formation of a paircd-H final states, a significant 
confirmation of aspects of some current influential mod- 
els^. Furthermore, we obtain final states with additional 
dangling bond defects, well separated from each other, 
again in agreement with key experiments^. Our Simula- 
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tions surfer from limitations: small cells and certainly a 
limited sampling of the possible motions of H, limited 
simulation times and approximations of various forms 
described below. Nevertheless, we believe that studies 
of this type offer significant promise as an "unbiased" 
means to discover the importance of H motion, its topo- 
logical and electronic consequences, and provide another 
needed piece to a complex puzzle. 

The rest of the paper is organized as follows. In Scc.lTTl 
we discuss the approximations used in the ab initio lo- 
cal basis code SIEST A 27 ? 28 ; 29 , describe the procedure 
for generating the models and also discuss the methods 
used to simulate both the electronic ground state and the 
light excited state. In Sec. lIIII we present detailed discus- 
sions of the molecular dynamics calculations of Hydro- 
gen dynamics and its consequences both in the electronic 
ground state and in the presence of light excitation. The 
change in the electronic structures and vibrational modes 
are explained in detail. We present conclusions in Sec. lIVI 



II. METHODOLOGY 

A. Total energies, electronic structure and 
dynamical simulation 

Simulation of complex a-Si:H models require accu- 
rate interatomic interactions (in particular, H energetics 
is highly delicate in a-Si:H) 30 - 31 . Therefore, our den- 
sity functional simulations were performed within the 
generalized gradient approximation 32 (GGA) or the lo- 
cal density approximation (LDA) using the first prin- 
ciples code SIESTASL 2 ^ 9 .. Norm conserving Troullier- 
MartinsiS pseudopotentials factorized in the Kleinman- 
Bylander— form were used. All calculations in this paper 
employed optimized double £ polarized basis sets (DZP), 
where two s and three p orbitals for the H valence elec- 
tron and two s, six p and five d orbitals for Si valence 
electrons were used. The structures were relaxed until 
the forces on each atom were less than 0.04 eV/A. We 
used a plane wave cutoff of 100 Ry for the grid (used 
for computing multi-center matrix elements) with 10~ 4 
for the tolerance of the density matrix in self consistency 
steps. We solved the self-consistent Kohn-Sham equa- 
tions by direct diagonalization of the Hamiltonian and a 
conventional mixing scheme. The T point was used to 
sample the Brillouin zone in all calculations. 

Density functional theory in the LDA, or with gradient 
corrections, maps the ground state many-electron prob- 
lem onto a system of non- interacting fermions 35 . In prin- 
ciple, the eigenvalues of the resulting single-particle equa- 



tions are not true excitation energies and the spectral gap 
between occupied and unoccupied states is well known 
to be incorrect. Nevertheless, the eigenvectors of the 
problem (the Kohn-Sham orbitals) have been shown to 
be very similar to quasiparticle states from "GW" calcu- 
lations, in which the self-energy is expressed as a product 
of the single particle Green's function "G" and the dy- 
namically screened Coulomb interaction "W" as used in 
many-body calculations 37 . For Si, C and LiCl, Hybertsen 
and Louie 3 ^ found 99.9% overlap between "GW" states 
and the Kohn-Sham orbitals. On an empirical level for 
amorphous materials, there are many indications that it 
is profitable to interpret the Kohn-Sham orbitals "liter- 
ally" for comparisons to experimen-fe 3 ^^. This provides 
some rationale for interpreting the Kohn-Sham orbitals 
as quasiparticle states, as we shall in some subsequent 
parts of this paper. 



B. Models 

A predictive simulation requires a physically plausi- 
ble model that represents the topology of the network 
and yields an accurate description for dynamics of the 
atoms. In this article we have used two different su- 
percell models: a 223 atom a-Si:H model (Model-I) and 
a 71 atom a-Si:H model (Model-II). These models are 
generated from a 64 atom and a 216 atom a-Si models 
which were generated by Barkema and Mousseau^i using 
an improved version of the Wooten, Winer, and Weaire 
(WWW) algorithm— respectively. 

To create the a-Si:H environment a) We started from 
a 216 atom a-Si model with two dangling bonds, we re- 
moved two silicon atoms resulting in the formation of 
additional vacancies. All of the vacancies except one are 
then terminated by placing a H atom at about 1.5 A from 
the corresponding Si atom. This yield a 223 atom Model- 
I. b) We started from a defect free 64 atom a-Si model, 
we removed three silicon atoms resulting in the forma- 
tion of vacancies. All of the dangling bonds are then 
terminated by placing a H atom at about 1.5 A from the 
corresponding Si atom to generate a 71 atom Model-II. 
We then repeated this supercell surgery at other sites to 
generate an ensemble of three configurations to obtain 
some insight into the formation of the structure and its 
bonding in solid state. Finally these newly generated 
structures are well relaxed using conjugate gradient op- 
timization technique. While such a procedure is clearly 
unphysical, it is worth pointing out that the resulting 
proton NMR second moments of the clusters created are 
similar to the broad component of the lineshape observed 
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in experiments . 



C. Excited state dynamics and promotion of 
carriers 

Defects in an amorphous network may lead to local- 
ized electron states in the optical gap or in the band 
tails. If such a system is exposed to band gap light, 
it becomes possible for the light to induce transitions 
from the occupied states to unoccupied states. For the 
present work we will not concern ourselves with the sub- 
tleties of how the EM field introduces the transition, we 
will simply assume that a photo-induced promotion oc- 
curs, by depleting the occupied states of one electron 
"forming a hole" and placing the electron near the bot- 
tom of the unoccupied "conduction" states. The idea is 
that a system initially at equilibrium will not be after 
the procedure: Hellmann-Feynman forces^ due to the 
occupation change will cause structural rearrangements, 
which may be negligible or dramatic, depending on the 
flexibility or stability of the network, and the localiza- 
tion of the states. The changes in force will initially be 
local to the region in which the orbitals are localized, fol- 
lowed by transport of the thermal energy. In general, it 
is necessary to investigate photo-structural changes aris- 
ing from various different initial and final states, though 
only well localized states near the gap have the poten- 
tial to induce structural change 22 - 45 ji£^£^&. The simu- 
lated light excited state is achieved by implementing: a) 
starting from the well relaxed model, we make the oc- 
cupation change by adding an additional electron just 
above the Fermi level, b) we keep the system in this 
excited state for lOps (20000 MD steps with time step 
r=0.5 fs between each MD steps), and maintain a con- 
stant temperature T=300K, c) after lOps, we put the 
system back into the ground state and relax to minimize 
the energy. The method has been described in additional 
detail elsewhere^. 



III. HYDROGEN DYNAMICS 

We have performed extensive MD simulations of net- 
work dynamics of a-Si:H both in an electronic ground 
state ("light-off") and a simulated light-excited states 
("light-on") for the two models, Model-I and Model-II 
described above. In the next sections, we present a de- 
tailed calculation of hydrogen diffusion, its mechanisms 
and consequences on the structural, electronic and vi- 
brational properties in both electronic ground state and 



light excited state. 

A. Hydrogen motion: Ground State 

To analyze the diffusion mechanism in the ground state 
we performed a MD simulation for five different temper- 
atures, and tracked the trajectories and bonding infor- 
mation of all the H and Si atoms in the network. In all 
the cases, the MD simulations show diffusion of hydro- 
gen in the cell and as a consequence, the network exhibits 
bond break and formation processes. The pattern of dif- 
fusion differs for individual H atoms depending upon the 
geometrical constraints around the diffusing H atom. 




FIG. 1: Trajectory for three different hydrogen atoms (H219 
and H220) in the ground state, which shows the diffusion and 
trapping of the atom for Model-I model. The total time for 
the trajectory is lOps. 

In order to characterize the trajectories of H diffusion 
in the ground state, we have selected two diffusive H 
atoms, (H219 and H220), and plotted their trajectories at 
T=300K in Fig. [TJ The trajectories for both H 2 i 9 and 
H220 atoms show diffusion in which the H atoms spend 
time being trapped in a small volume of the cell which 
is followed by rapid emission to another trapping site. 
In order to examine how the bond rearrangement takes 
place in the network while the H atom is diffusing, we 
tracked each hydrogen atoms and computed its bonding 
statistics. 

In Fig. |21 we show the Si-H bond length between one 
of the diffusing H atoms (namely H 2 ig) and relevant Si 
atoms (Sigo and Sii2s) with which it forms a bond while 
diffusing and Si2os ■ As we can see from Fig. in the 




FIG. 2: The Si-H bond length between the diffusing H (H219) 
and three different Si atoms, (Si9o, S1208, and Sii28 as a func- 
tion of time in the electronic ground state for Model-I. The 
total time for the trajectory is lOps. 

first 4ps H219 is bonded with Sigo with a bond length of 
1.5 A and trapped for a while until it breaks and hops 
to form another bond with Sii28- In the first ~4 ps, the 
bond length between H219 and Sii28 fluctuates between 
3.8 A and 2.5 A. However, after ~4 ps we observed a 
swift bond changes in a very short period of time ~0.1 
ps when the H219 atom comes out of the trapping site 
and hops to form a bond with Sii28 and trapped there 
for a very long time period of ^6 ps. This process of 
trapping and hopping is typical for the highly diffusive H 
atoms. 

To study atomic diffusion we computed the time aver- 
age mean squared displacement for both H and Si atoms 
for a given temperature using 

(a 2 (a, T)) Ume = j^j^ E E ~ f * Q (°)| 2 ' 

t—1 i—1 

(!) 

where the sum is over particular atomic species a (Si or 
H), N a and fi a (t) are total number and coordinates of 
the atomic species a at time t respectively, and Nmd is 
the total number of MD steps. 

The time average mean square displacement for Model- 
II for five different temperatures was calculated using 
Eq. Q for H atoms in the supercell in the electronic 
ground state (light-off) and it is shown in Fig. [3J We 



FIG. 3: Time average mean square displacement for H as 
a function of temperature of MD simulation in electronic 
ground state for Model-II model . 

have observed a strong temperature dependence of H dif- 
fusion. This result will help us to compare the diffusion 
of H in the electronic ground state with the light excited 
state to be discussed in the next section. 



B. Hydrogen diffusion: light excited state 

Similar to the case of electronic ground state, we an- 
alyzed the diffusion of H in the light excited state by 
performing a MD simulation. We tracked the trajec- 
tories and bonding statistics of Si and H atoms in the 
supercell. Our MD simulation in the light excited state 
show enhanced hydrogen diffusion and consequently in- 
creased bond breaking and formation that leads to struc- 
tural changes in the network. 

For the purpose of analyzing the difference in the dif- 
fusion mechanism of H in the light excited state case 
as compared with the ground state, we performed simi- 
lar calculations described in the previous sections for the 
light excited state case. To see the trajectories of H in 
the light excited state, we have again selected two diffu- 
sive H atoms, (H219 and H220) from the larger Model-I, 
and plotted their trajectories in the light excited state 
in Fig. 21 The trajectories show the diffusion of H in 
the presence of different trapping centers, a region where 
the H atom spends more time before it hops and moves 
to another trapping site. However, in this case we ob- 





FIG. 4: Trajectory for three different hydrogen atoms (H219 
and H220) which shows the diffusion and trapping of the atom 
for Model-I in the light excited state. The total time for the 
trajectory is lOps. 



served enhanced diffusion and more trapping sites and 
hopping. These trapping and hopping processes continue 
until two hydrogens form a bond to a single Si atom to 
form a metastable S1H2 conformation or until two hydro- 
gens form a bond to (a) two different Si atoms which are 
bonded to each other, to form (H-Si-Si-H) structure or 
(b) two different Si atoms which are not bonded but close 
to each other to form (H-Si Si-H) structure. This is in 
agreement with a basic event of the H collision model 6 
and other H-pairing models 21 . 

By tracking each H atom, we computed its bond- 
ing statistics and examine the bond rearrangements. In 
Fig. [3] we show Si-H bond length as a function of time 
between one of the diffusing H atoms (H219) and three 
other Si atoms (Sigo, Sii28> an d Si2os) with which it forms 
a bond while diffusing in the network. The initial trap- 
ping time, where H219 is bonded with Sigo, is reduced to 
~1.8 ps when the light is on from ~4 ps when the light is 
off. This is followed by another trapping site where H 2 ig 
is bonded with Si2os for another ~2.1 ps. The H 2 ig hops 
out of the trapping site and forms a bond with Sii28 and 
trapped for ~4.3 ps before it finally hops out from the 
trapping site and forms another bond with Si2os where 
it gets trapped again and form a silicon dihydride (SiH.2) 
structure. As we can see from Fig- El the pattern of diffu- 
sion is quite different from the ground state: In the light 
excited state case we observed a) more number of trap- 
ping sites and less trapping time with frequent hopping, 
b) enhanced hydrogen diffusion, and c) increasing num- 
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FIG. 5: The Si-H bond length between the diffusing H (H219) 
and three different Si atoms (Sigo, Siias, and Si2os) with which 
H219 forms a bond (one at a time) while it is diffusing as a 
function of time for Model-I, in the light excited state. The 
total time for the trajectory is lOps. . 



ber of bond rearrangements and newly formed dihydride 
structural units. 

The atomic diffusion in the light excited state case has 
also been examined using the time average mean squared 
displacement for both H and Si atoms for different tem- 
peratures using Eq. ^ for both Model-I and Model-II. 
The results from Model-II is shown in Fig. HJ] For all 
the temperatures considered, our simulation results show 
enhanced diffusion of Hydrogen for the case when the 
light is "on" as compared with the case where the light 
is "off" . Consistent with the work of Isoya- 6 -, the hop- 
ping of H is apparently stimulated by the electron-hole 
pair. The enhanced diffusive motion of H in the photo 
excited state relative to the electronic ground state arises 
from the strong electron-lattice interaction of the amor- 
phous network, and an effect of "local heating" and sub- 
sequent thermal diffusion^ 6 ' initially in the spatial volume 
in which the state is localized. The same calculations 
has also been performed on the larger model Model-I at 

T=300K in which, the time average mean square dis- 
2 2 

placement for H is 2.66 A for the light excited state and 
1.10 A for the electronic ground state. These results 
again show and confirm an enhanced hydrogen diffusion 
for the case of light excited state. In all the cases no 
enhanced motion for Si is observed 



6 



H (light-on) 




8* 0.6 - 



^ 0.4- 
0.2 - 



200 400 600 800 

Temperature (K) 

FIG. 6: Time average mean square displacement for Has a 
function of temperature of MD simulation in the light excited 
for Model-II . 



C. Consequences of Hydrogen diffusion 

1. Formation of dihydride structure 

In the two scenarios that we considered, MD simula- 
tion in electronic ground state (light off) and simulated 
light-excited state (light on) we have observed an impor- 
tant difference. In the light-excited state, in addition to 
bond rearrangements and enhanced hydrogen diffusion, 
we have observed a preferential formation of new struc- 
ture: S1H2, with an average distance of 2.39 A for the 
pair of hydrogens in the structure, (H-Si-Si-H) and (H-Si 
Si-H) with H-H separation which ranges from 1.8 A to 
4.5 A. However, in the electronic ground state, we have 
obtained rearrangement of atoms including hydrogen dif- 
fusion, without formation of SiH2 structure in the super- 
cell. The mechanisms for the formation of these struc- 
tures in the light-excited state follows breaking of H atom 
from Si-H bond close to the dangling bonds and diffusion 
to the nearest weakly bonded interstitial sites (or dan- 
gling bonds). This mobile H atom then collides (forms 
a metastable bond) with another Si+DB structure or 
breaks an Si-Si bond to form another Si-H bond. This is 
attributed to the fact that the dangling bond site is mov- 
ing to accommodate the change in force caused by the 
additional carrier and also because hydrogen is moving 
through weakly bonded interstitial sites with low activa- 



tion barrier for diffusion until it is trapped by a defect*^ 
We find that there are two different modes of bond for- 
mation for the mobile hydrogen. The first is when two 
mobile hydrogen atoms, H m , collide with two Si atoms 
and form a metastable (H-Si-Si-H) or (H-Si Si-H) struc- 
ture and the second one is when the mobile hydrogen 
moves until it encounters a preexisting Si-H+DB struc- 
ture and makes a bond to form a SiH2 structure. 

Consequently, our calculations show two basic ideas 
for the diffusion of H in the light-excited state: 1) the 
diffusion of hydrogen doesn't only break a Si-H bond but 
it also breaks a Si-Si bond and 2) the possibility that two 
mobile H atoms might form a bond to a single Si atom 
to form a metastable SiH2 structure in addition to the 
formation of (H-Si-Si-H) and (H-Si Si-H) structures. 

In the Model-II, the two hydrogens involved in the for- 
mation of the S1H2 structure initially were 5.50 A apart 
and bonded to two different Si atoms (Si-H) which were 
separated by 4.86 A. With thermal simulation in the light 
excited state, the two hydrogen atoms dissociate from 
their original Si atoms and becomes mobile until they 
form the SiH2 structure, in which the H-H distance be- 
comes 2.39 A. We have observed similar pattern of H 
diffusion, bond rearrangements and formation of S1H2 
structure near the DB for the other two configurations 
considered in the simulation. The same phenomenon is 
observed in the case of Model-I. The two hydrogens in- 
volved in the formation of the SiH2 structure initially 
were 3.29 A apart and bonded to two different Si atoms 
(Si-H) which were separated by 3.92 A. With thermal 
simulation in the light excited state, the two hydrogen 
atoms dissociate from their original host and becomes 
mobile until they form the SiH2 structure, in which the 
H-H distance becomes 2.45 A. We have summarized the 
results that show before and after MD calculations of 
H-H distance (in SiH2 structure) for Model-I and three 
different configurations of Model-II in the case of light 
excited state in Table H] 

2. Change in the electronic properties 

In order to understand the electron localization we 
used the inverse participation ratio, 2", 

N 

l = Y / h(E)} 2 (2) 

i=l 

where qi{E) is the Mulliken charge residing at an atomic 
site i for an eigenstate with eigenvalue E that satisfies 
Yh )] = 1 ano - N i s t ne total number of atoms in the 
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TABLE I: The H-H distance in the SiH2 configurations and 
the Fermi energy of the system before and after MD simula- 
tions in the light excited case. 



Config- 
-urations 



H-H distance 
before MD after MD 

(A) (A) 



1 (Model-II) 

2 (Model-II) 

3 (Model-II) 

4 (Model-I) 
Average 



5.50 
3.79 
4.52 
3.29 



2.39 
2.36 
2.36 
2.45 
2.39 



cell. For an ideally localized state, only one atomic site 
contributes all the charge and so I = 1. For a uniformly 
extended state, the Mulliken charge contribution per site 
is uniform and equals 1/N and so X — 1/N. Thus, large 
X corresponds to localized states. With this measure, 
we observe a highly localized state near and below the 
Fermi level and a less localized state near and above the 
Fermi level. These states, highest occupied molecular or- 
bitals (HOMO) and lowest unoccupied molecular orbitals 
(LUMO), are centered at the two dangling bonds in the 
initial configuration of the model. The energy splitting 
between the HOMO and LUMO states is 1.08 eV. Figure 
[7| (a) shows the Fermi level and X of these two states and 
other states as a function of energy eigenvalues in the 
relaxed electronic ground state. 

This picture changes when we excite the system and 
perform a MD calculation in which we observe enhanced 
diffusion of hydrogen and subsequent breaking and for- 
mation of bonds. Since electron-phonon coupling is large 
for localized states^, the change of occupation causes 
the forces in the localization volume associated with the 
DB to change and the system moves to accommodate the 
changed force. Consequently, the hydrogen atoms close 
to the DB sites start to move in the vicinity of these 
defects either to terminate the old DB's or to break a 
weak Si-Si bond and by doing so, create new DB de- 
fects on nearby sites. As shown in Fig. (b) we observe 
the formation a highly localized state and appearance of 
three less localized states, that correspond to the newly 
formed defect levels after simulated light-soaking. These 
processes induce transition of electrons from the top of 
the occupied states to the low-lying unoccupied states 
which is reflected in the smaller value of X for the initial 
HOMO state and an increase in the X for the LUMO 
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FIG. 7: (Color Online) The inverse participation ratio I of 
the eigenstates versus the energy eigenvalues, (a) in the re- 
laxed electronic ground state and (b) in the relaxed simulated 
light-excited state (light excited MD followed by relaxation), 
with their respective Fermi energy in the first configuration 
of relaxed Model-II. The inset (c) shows the electron density 
of states with the Fermi level shifted to zero for the relaxed 
simulated light-excited state. 



state. 

The X of the HOMO state, where the state is ini- 
tially localized, decreases from 0.158 to 0.060 after photo- 
excitation, while the X of the LUMO state increases from 
0.045 to 0.142. The splitting energy between the HOMO 
and LUMO states has also declined to 0.723 eV. The 
newly formed defects with lower energy splitting between 
the HOMO and LUMO states suggest a presence of car- 
rier induced bond rearrangements in the supercell. The 
comparisons for the energy and X of the system before 
MD (as relaxed) and after MD is given in Table ILTl 

In addition, analysis of the spatial distribution of the 
configurations shows that the H atoms close to the dan- 
gling bonds (< 4.0 A) are most diffusive and the Si atoms 
which make most of the bond rearrangements includ- 
ing the Si atom in the SiH2 configurations are close (< 
5.50 A) to the dangling bonds. These show the additional 
charge carrier induces change in the forces around the 
dangling bonds and consequently rearranges the atoms 



s 



60 



(a) 




0.04 



0.03 



0.02 



0.01 



0.02 



0.01 



-10 12 
Energy (eV) 

FIG. 8: (Color Online) The energy density of states and the 
inverse participation ratio X of the eigenstates versus the en- 
ergy eigenvalues, (a) in the relaxed electronic ground state 
and (b) in the relaxed simulated light-excited state (light ex- 
cited MD followed by relaxation), both the electron density 
of states and the inverse participation ration are plotted with 
the Fermi level shifted to zero. 



influence these numbers, but it is clear that the defect 
creation is not very local because of the high diffusivity 
of the H. 

The same calculation have been performed on Model- 
II. In Fig. Owe have plotted both energy density of states 
and inverse participation ratio as a function of energy in 
the light excited state case before and after the MD sim- 
ulation. As can be seen from the figure we obtained more 
localized states in the middle of the gap which are caused 
due to an increase in the number of defects upon light 
excitation. This supports that the diffusion of hydro- 
gen not only forms preferential dihydride structures but 
also increase the number of defects in agreement with our 
findings for the smaller cell model of asiH-71. 



3. Change in the vibrational properties 

For an amorphous solid, the vibrational density of state 
is a sum of 3iV (N is the number of atoms) delta functions 
corresponding to the allowed frequency modes. Starting 
with the relaxed Model-II subsequent to MD in the light 
excited state, we computed the vibrational energies (vi- 
brational modes) from the dynamical matrix, which is 
determined by displacing each atom by 0.02 A in three 
orthogonal directions and then performing ab initio force 
calculations for all the atoms for each displacement to ob- 
tain the force constant matrix, and with diagonalization, 
phonon frequencies and modes. 



TABLE II: The energy and the inverse participation ratio 1 
of localized states HOMO, LUMO, LUMO+1 and LUMO+2 
before and after the MD for Model-II. 



Eigenvalue T 



TABLE III: Frequency for some of the Si-H vibrational modes 
of the SiH2 conformation for the first two configurations of 
the Model-II obtained from our MD simulations and their 
corresponding experimental values-"*—*—. 



before MD after MD before MD after MD 





(eV) 


(eV) 






HOMO 


-4.32 


-4.40 


0.158 


0.060 


LUMO 


-3.24 


-3.68 


0.045 


0.142 


LUMO+1 


-2.88 


-3.08 


0.037 


0.064 


LUMO+2 


-2.66 


-2.87 


0.030 


0.042 





Rocking 


Scissors 


Stretch 


Configurations 


cm -1 


cm -1 


cm -1 


1 


629 


810 


2025 


2 


625 


706 


2047 


Experiment 


630 


875 


2090 



around the dangling bond sites and eventually forming 
an SiH2 structure. On average the newly formed defect 
sites are 3.80 A and 4.70 A far away from the two ini- 
tial defect sites. The newly formed SiH2 structure is (on 
average 4.11 A away from the initial defect sites. It is 
probable that limitations in both length and time scales 



In our calculations, the VDOS shows H modes of vi- 
brations in the range (600-900) cm -1 and also in the 
range (1800-2100) cm -1 . We have examined the vibra- 
tional modes to pick out those modes arising only from 
SiH2. We reproduce the vibrational modes of SiH2 and 
their corresponding experimental values^i^S in Table 
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IIIII The first mode is the rocking mode at 629 cm -1 and 
625 cm -1 ; the second is the scissors mode at 810 cm -1 
and 706 cm -1 and the last is the asymmetric stretching 
mode that occur at 2025 cm -1 and 2047 cm -1 for the 
first and second configurations respectively. These re- 
sults are in good agreement with the IR absorption spec- 
tra for the SiFf2 structure. The comparison of our results 
for the vibrational modes of SiH2 with the experiment is 
summarized in Tabic ITTT1 The results shown in Table ITTT1 
are sensitive to the basis sets used in the calculation, in 
agreement with other work emphasizing the delicacy of 
H dynamics^. 

IV. CONCLUSION 

We have presented a direct ab-initio calculation of 
network dynamics and diffusion both for the electronic 
ground state and light-excited state for a-Si:H. We com- 
puted the preferential diffusion pathways of hydrogen 
in the presence of photo-excited carriers. In the light- 
excited state, we observe enhanced hydrogen diffusion 
and formation of new silicon dihydride configurations, 
(H-Si-Si-H), (H-Si Si-H), and SiH 2 . The two hydrogens 
in the SiH2 unit show an average proton separation of 



2.39 A. The results are consistent (a) with the recent 
NMR experiments and our previous studies, and (b) with 
the hydrogen collision model of Branz and other paired 
hydrogen model in the basic diffusion mechanism and for- 
mation of dihydride structures. In contrast, simulations 
in the electronic ground state do not exhibit the ten- 
dency to SiH2 formation. Undoubtedly, other H diffusion 
pathways exist, and the importance of larger simulation 
length and time scales as well as effects of promotions 
involving different states (which could include strain de- 
fects and floating bonds^) should be undertaken. For the 
first time, we show the detailed dynamic pathways that 
arise from light-induced occupation changes, and provide 
one explicit example of defect creation and paired H for- 
mation. 
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